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Abstract 


Two finite element based computational methods, Smoothing Element 
Analysis (SEA) and the inverse Finite Element Method (iFEM), are 
reviewed, and examples of their use for structural health monitoring are 
discussed. Due to their versatility, robustness, and computational 
efficiency, the methods are well suited for real-time structural health 
monitoring of future space vehicles, large space structures, and habitats. 
The methods may be effectively employed to enable real-time processing 
of sensing information, specifically for identifying three-dimensional 
deformed structural shapes as well as the internal loads. In addition, 
they may be used in conjunction with evolutionary algorithms to design 
optimally distributed sensors. These computational tools have 
demonstrated substantial promise for utilization in future Structural 
Health Management (SHM) systems. 


Introduction 

The Integrated Vehicle Health Management (IVHM) research at NASA Langley Research Center 
explores safe, reliable, and affordable technologies for NASA’s future long duration space missions [1]. 
Integral to this research is mitigation of aerospace vehicle accidents due to structural failures. For long 
duration space missions, real-time monitoring of structural, propulsion, thermal protection, and other 
critical systems will be required. To achieve such capabilities, space vehicles and habitation structures 
will be designed with diverse arrays of optimally distributed in-situ sensors. The sensing technologies will 
be part of advanced data systems architectures that will process, communicate, and store massive amounts 
of SHM data. Special-purpose structural analysis and design algorithms will be necessary to incorporate 
SHM sensing data for the diagnosis and prognosis of structural integrity, and for the purpose of optimal 
design and construction of such structures. Structural integrity information will be utilized within an 
IVHM system, resulting in the safe and effective operational vehicle management and mission control. 

In recent years, ample focus has been on advancing the state-of-the-art in sensing technologies and signal 
analysis. These capabilities are designed to acquire accurate structural response information and to infer 
the state of structural deformation and potential damage and defects. Much of this research has been in the 
area of fiber Bragg grating (FBG) optical sensors. When bonded to or imbedded in load-carrying 
structures, FBG sensors may provide high-quality multi-point strain measurements in a variety of 
applications including civil, marine, and aerospace. A key step in analyzing strain data is to infer or 
reconstruct an accurate representation of the deformed structural shape. Mathematically, this type of 
analysis represents an ill-posed inverse problem and has been referred to as shape sensing. Several shape- 
sensing algorithms have been explored and thus far applied to relatively simple structures [2-5]. FBG 
optical sensors provide lightweight distributed capabilities for performing shape sensing computations 
which are essential in facilitating digital control of aerodynamic surfaces during flight. This is particularly 
relevant to flexible-wing vehicles, such as a Helios class of aircraft shown in Figure 1, requiring 
automated procedures to control wing dihedral in flight [6]. 




Figure 1. Helios solar powered aircraft during first flight on July 14, 2001 at Dryden Flight Research Center, 
California. 

In general, the class of Unmanned Aerial Vehicles (UAV) may derive substantial performance benefits 
using real-time wing surface control systems. For large space structures, including solar sails and 
membrane antennas, knowing the current three-dimensional shape of the structure may maximize 
spacecraft performance [7-8]. 

This paper reviews two notable full-field computational algorithms, known as SEA [9-10] and iFEM [4- 
5]. The methods combine the attributes of computational efficiency, versatility, and robustness that are 
necessary for real-time, large-scale SHM applications. The algorithms may serve as both design tools for 
the development of optimal sensing systems and enabling tools for the processing of sensing information 
in real time. Following the mathematical description of SEA and iFEM and discussions of their salient 
features, several SHM applications of the two methods are highlighted. These include (1) shape sensing 
and structural anomaly detection in beam structures instrumented with FBG sensors, (2) identification of 
delamination damage in composite laminates, and (3) design of optimal strain sensor locations. 

Smoothing Element Analysis 

Smoothing Element Analysis (SEA) provides a robust and versatile computational capability for 
smoothing discrete data of either experimental or numerical origin that are defined over arbitrary 
geometric domains. In addition to obtaining smooth fields, the method computes smooth first-order 
derivatives of the field of interest. SEA is a finite element method derived on the basis of a Penalized 
Discrete Least Squares (PDLS) variational principle. The method treats any discrete data as 
corresponding to a scalar quantity. 

The general problem is that of smoothing a set of experimental, inherently noisy measurements obtained 
by an array of sensors that measure a certain response quantity. The structure may possess complex 
geometric features and may undergo general deformations (linear or nonlinear) under the action of 
applied loads and imposed boundary restraints. In Figure 2(a), the applied loads, boundary restraints, and 
sensor locations are identified on a structural model. A corresponding candidate SEA mesh for the 
instrumented structure is shown in Figure 2(b). The SEA mesh topology is strongly dependent upon the 
number of sensors, their locations, and the element polynomial interpolations of the field variables. 

Consider a planar surface where numerous sensors measure the response quantity of interest. The surface 
is confined to a Cartesian frame of reference and is discretized with a mesh consisting of ‘smoothing’ 
elements, as shown in Figure 2. For a single smoothing element defined over the area A e that 

encompasses n points of measured data (for example, strain 8^ , where q = l,...,n), the PDLS error 
functional may be expressed as 
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where e(x q ) a smooth strain is interpolated over the element and evaluated at the sensor position 

x q = (x q ,y q ); 0 x and 0 y are independent variables that, respectively, approach very closely the partial 

derivatives 8 and e (refer to Eqs. (2)). In general, the symbol 8 may represent an arbitrary quantity of 

interest, e.g., displacement or temperature. For superior numerical accuracy, the complete array of 
measured data must be normalized with respect to its maximum magnitude. 

The first term in Eq. (1) is weighted least squares with the weighting coefficient defined as the inverse of 
the number of sensors (obviously, other normalization options exist). Thus, the measured data are 
weighted equally within an element. The second term in Eq. (1) is a penalty constraint functional which, 
for a relatively large value of the dimensionless parameter a (~10 6 ), ensures the limiting condition of C 1 
continuity for 8, 

£ ,x^ 0 x and ^,y * 0 y in A e (2) 

where it is assumed that 8,0 x and 0 y are interpolated using C°-continuous shape functions. The third 
term in Eq. (1) is a regularization function in which (3 is a positive regularization parameter. Assigning 
small values of (3 (~10 6 ) provides additional constraint conditions for the 0 x and 0 y variables. This 

yields nonsingular least squares equations even when some elements do not encompass any sensors or 
data. 

In accordance with the standard finite element procedures, the total PDLS error functional is computed as 
the sum of all element contributions and is minimized with respect to the nodal degrees-of-freedom 
associated with the 8, 0 x , and 0 y variables. The resulting linear equations are factorized, yielding point- 
wise smooth solutions for these variables. 



Figure 2. (a) A schematic of a structural model showing the load, boundary restraints, and sensor locations, (b) SEA 
mesh discretizing the structural model and encompassing the sensor positions. 
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The practical aspects of applying SEA include: (1) boundary conditions may be imposed on both the 
smoothed quantity 8 and the variables, 0 X , and 0 y , thus taking advantage of certain physical conditions 

that are known a priori, (2) higher-order derivatives of the measured quantity can be computed by 
applying the method sequentially in multiple steps, and (3) multiple sets of data may be smoothed 
efficiently, especially when the data are measured at the same sensor locations (e.g., three strain 
components measured by strain rosettes). For such cases, an efficient factorization algorithm of linear 
equations with multiple right-hand sides is required. 

Inverse Finite Element Method 

The iFEM has been developed for plate and shell structures [4-5] and is readily reduced to beams and 
frames. The method is formulated for the purposes of (1) reconstructing a smooth displacement field from 
strain data measured by in-situ strain sensors, and (2) computing the full-field strain, stress, and internal 
load (stress resultant) fields using the reconstructed displacement field. This part of the analysis follows 
the same basic procedures of the forward (direct) FEM. 

The first (and primary) part of iFEM performs a smoothing operation of all components of measured 
strain simultaneously. This requires all strain-tensor relations pertinent to the analytic theory be 
incorporated in the formulation. The iFEM formulation is based upon minimization of a least squares 
functional that uses Mindlin theory assumptions. The error functional for an inverse shell element of area 
A e that encompasses n strain-sensor locations may be stated as 

<&«(u) = ||e(u) -e e || 2 + ||k(u) -k e f + A ||g(u) -g e f 

where the squared norms are least-squares difference terms defined as 

||e ( u ) - e' If - T [ e (“), - ej ] 2 d\ dy 

H u )- k 1l ! -^f<2[ k < u ><- k '] dx dy (4) 

^ e i=\ 

||g (u ) - g £ || 2 - fy [g (u),. - g- ] 2 dx dy 

and A (~10 -6 ) denotes a dimensionless regularization parameter. The strain measures are defined in a 
customary manner in terms of C°-continuous kinematic variables u T = ju, v, w, 0 X , 0 y j as 

e( u ) T = {u x , v y , 0 X y + 0 y x } (In-plane strains) 

k (u) T = {0 y x , 0 x y , 0 X>X + 0 y y } (Curvatures) (5) 

g( u ) T = { w x + 0 y > w y + 0 X } (Transverse shear strains) 


4 



where u = u(x,y) and v = v(x,y) are the mid-plane displacements in the x and y directions, 
respectively; 0 x = 0 x (x, y) and 0 y = 0 y (x, y) are the rotations of the normal about the negative x and 
positive y axes, respectively; and w = w(x, y) is the deflection variable which is constant across the 
thickness coordinate z G [— t, t] , with 2t denoting the shell thickness. 


The e £ and k £ strain measures are readily derived from the measured surface strains on the top (+) and 
bottom (-) surfaces (refer to Figure 3) 
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The transverse shear strains g £ cannot be directly obtained from surface measurements. For thin shells, 
however, the g £ strains are negligibly small and therefore can be readily omitted from the formulation. In 
this case, the third term in Eq. (4) may be stated as the L 2 norm 

||g(u)-g e || 2 = / A g(u) 2 dxdy 

e \') 

The above term will therefore enforce, in a limiting sense, the Kirchhoff condition of zero transverse 
shear strains, 

g(u) T ={w >x + e y ,w >y + e x }-*o (8) 


Thus, when strains are measured on the top and bottom surfaces (or any other two surfaces parallel to the 
middle surface), the Kirchhoff-type (classical) shell behavior involving both stretching and bending may 
be modeled in a straightforward manner. For most practical applications this is sufficient to reconstruct 
accurately the deformation field of a given structure without the need to perform additional and costly 
computations to estimate the shear strain 

Because all strain-displacement relations are fulfilled, the strain compatibility relations of Mindlin theory 
are explicitly satisfied. Thus, given the distribution of strain sensors and strain values, the formulation 
provides a variationally consistent means of integrating the strains to obtain a smooth displacement field. 
There are no material properties being used to arrive at the displacement field since none exist in the 
strain-displacement relations. 



Figure 3. A panel instrumented with strain sensors. 
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As in the direct FEM, the error functional of Eq. (3) provides a mathematical framework for the 
development of a number of different finite element approximations. One such element, which is simple, 
computationally efficient, and robust, is a three-node, inverse shell element developed in [2]. The element 
has the typical shell kinematic description with six engineering degrees-of-freedom (dof s) defined at 
each node (see the illustration in Figure 4). To enable robust modeling of built-up shell structures, a 
drilling rotation degree-of-freedom, 0 z , is added at the element level. 

Following a suitable element interpolation scheme for the kinematic variables in the u vector, the 
standard finite element procedure is followed: (1) minimization of Eq. (3) with respect to the 
displacement dof s, (2) assembly of the element matrices into the global system of equations, using the 
appropriate element-to-global displacement transformations, and (3) imposition of the kinematic 
boundary conditions consistent with the physical constraints on the u displacements. 

The resulting global system of algebraic equations has the form 

A d = f (9) 

where d is a vector of all active kinematic dof s in the discretization; A = A(Xj) is a square matrix that, 
for a given element definition, mesh and structural topology, is only a function of the strain-sensor 
positions, x i . Provided that at least a single strain sensor resides within an element and the appropriate 
kinematic boundary conditions are imposed, the A matrix is nonsingular. The right-hand side vector 
f = f(E + ,S-) is a function of the strain values only. Factorization of Eq. (9) results in a unique solution 
for the d displacement vector. Subsequently, C° continuous element displacements u are fully 
determined, thus ensuring C u displacement continuity across the entire domain of the discretized 
structure. 



y>v 



Figure 4. A three-node inverse shell element. 
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Full-Field Strains and Internal Loads 


Once the element displacement field u is obtained, the element level ‘recovery’ of strains, stresses and 
stress resultants follows standard finite element procedures. First, the in-plane and curvature strain 
measures are computed from Eqs. (5), leading to the computation of the strains 

Fxx> Sy’ Yxy} T s e(u) + zk(u) (10) 


Assuming linear elastic behavior, the stresses O = |a xx , cr yy , T xy j are obtained from the stress-strain 
relations 


a = D{e xx , £ yy , Y xy } T (11) 

The stress resultants R = |N X? N y , N xy , M x? M y , M xy | are computed from the constitutive relations 
that involve the strain measures (see Figure 3) 



where D and Q are the appropriate constitutive matrices. The transverse shear forces can then be 
calculated from the equilibrium equations of plate theory, 


V. =M. 


+ M xy , y and V y = M xy 


+ M„ 


y>y 


(13) 


Computational Speed and Real-Time Applications 


The issue of computational speed is paramount for real-time applications. An on-board computer should 
be able to perform the computational tasks fast enough to provide the requisite response data to an 
automated control system and the pilot visualization display. The special features of iFEM appear 
perfectly suited for such an environment. The iFEM computational paradigm is as follows. The iFEM 
structural model of an aerospace structure is fully defined in terms of its geometry, discretization by the 
inverse finite elements, the displacement boundary conditions, and the strain sensors positions. Since the 
system matrix A is a function of the strain-sensor positions only, the matrix is inverted once and for all, 
and remains unchanged during the operation of the vehicle. Thus, only the right-hand side vector f is 
recomputed in real time (since the measured strains change during flight) followed by the fast dot-product 
operation 


d = A f 


(14) 


The displacement data d (i.e., the deformed shape) is then fed into an automated control system to 
perform the necessary control-surface adjustments consistent with the current deformed shape of the 
structure, e.g., a wing. Furthermore, the recovered internal loads, computed according to Eqs. (11)-(13), 


7 



may be displayed on the pilot’s monitor to ascertain the current internal loads environment on the vehicle. 
Using the internal loads information, appropriate failure criteria may be calculated, stored, and displayed 
in real-time. 

Meshing Strategies and Optimization of Sensor Locations 

Both SEA and iFEM computational tools offer various meshing possibilities to be explored for a given 
distribution of sensors. In Figure 5, two basic meshing strategies are demonstrated on an idealized wing 
model with an arbitrary distribution of strain sensors. The first mesh is labeled 1-to-l, meaning that a 
single sensor resides within one element. The second mesh is labeled n-to-1 signifying that n sensors are 
positioned within an element, where n is an arbitrary number for each element. 

The three-node iFEM shell element previously described has constant in-plane and bending strains. Thus, 
to construct the highest-density mesh, that would be optimal from the accuracy standpoint, a single sensor 
would need to be positioned at the centroid of each element. In general, the three-node element 
performance is robust enough so that no specific restriction on the location of the sensors is required. 

Application of the SEA and iFEM tools may also facilitate the design of optimal distribution of sensors 
on the structure. Herein, for a fixed number of sensors, the aim is to obtain a particular distribution of 
sensors that would provide the most accurate response representation reproduced by a given 
computational tool at hand. Recently, this issue was explored by way of numerical studies of statically 
loaded plate structures using an iterative procedure involving Genetic Algorithms (GA), FEM, and SEA 
[11]. A procedure for utilizing the iFEM algorithm within the GA optimization scheme was also 
examined. 





Figure 5. Idealized wing model with an arbitrary distribution of strain sensors, and two alternative meshing 
strategies with respect to the locations of sensors. 

Shape Sensing Analysis using iFEM for Built-up Shell Structures 

In Figure 6, graphical results for a computational validation of iFEM are shown for a typical aerospace 
built-up structural component. The hat-stiffened panel with a circular cutout and fully clamped along the 
left and right edges is subjected to a concentrated force as shown. The ‘experimental’ strain 
measurements are obtained by means of the direct FEM (ABAQUS) shell solution with the strains 


computed at the element centroids. These strains are then seeded into an iFEM model having the same 
mesh as the FEM model. The resulting displacement field produced by the iFEM solution is less than 
0.1% in error compared to the FEM displacements; hence the displacement results are graphically 
indistinguishable (see the detail in Figure 6). Moreover, provided both the FEM and iFEM models belong 
to the same approximation space (i.e., they are based on the same underlying structural theory, mesh 
refinement, and element interpolations), the iFEM model would reconstruct the FEM fields 
(displacements, strains and stresses) exactly. 

Hat-Stiffened Panel 



Figure 6. Hat-stiffened panel deformed under a concentrated force with a comparison of FEM shell analysis and 
iFEM deformed shape reconstruction (‘shape sensing’) analysis. 

Shape Sensing and Structural Anomaly Detection using Fiber Optic Strain 
Systems 

The principles of Optical Frequency-Domain Reflectometry (OFDR) have recently been implemented at 
NASA Langley Research Center to provide a large-scale, dense distribution of strain sensors using optical 
fiber embedded with Bragg gratings. The theory and development of this technology, called the Fiber 
Optic Strain System (FOSS) are described in [12]. The algorithms and methods for deriving strain using 
OFDR measurements are documented in [13]. The fiber optic Bragg grating technology enables the 
distribution of thousands of sensors that are immune to moisture and electromagnetic interference and 
have negligible weight penalty. This technology can be conveniently coupled with the SEA- and iFEM- 
based computations that do not require strain measurements on a regular grid and provide robust 
interpolation and smoothing of the data relative to arbitrary locations (see the schematic in Figure 7). 



Figure 7. FOSS system coupled with SEA/ iFEM full-field reconstruction software. 
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The FOSS and iFEM technologies were successfully applied to a series of laboratory tests conducted on 
an aluminum bar undergoing bending deformations [14,15]. The first study focused on recovering the 
correct deformed shape (i.e., shape sensing analysis) of a beam based upon the discrete strain 
measurements provided by the FOSS system. In the second study, the aim was to detect structural 
anomalies represented by open holes. The objective was to identify structural anomalies that result in 
observable changes in localized strain but do not impact the overall surface shape. The iFEM analyses 
provided the full-field displacements and internal loads using strain data from in-situ fiber-optic sensors. 
Issues related to improvements in sensor application techniques and optimal levels of sensor density and 
layout schemes are yet to be investigated. 

Identification of Delamination Damage in Composite Plates 

In a recent computational study, a non-destructive detection of delamination damage in composite plates 
was addressed using high-frequency dynamic excitations [16]. The structural response was modeled by 
the finite element method using Mindlin-type plate elements. The approach does not require vibration 
measurements on the undamaged structure. The principal parameter used in assessing the state of 
delamination damage is the curvature tensor of the damaged structure. Once properly analyzed, the 
complete curvature tensor provides useful information for identifying the location and size of 
delaminations. 

The smoothed curvature data, corresponding to the high-frequency modes, serve as a measure of 
comparison for the raw mode-shape data that are obtained for a damaged structure. Whereas the modal 
shapes of the undamaged plate are presumed to be unavailable, the smoothed curvatures obtained via SEA 
provide the proper reference curvature fields. The implementation of a damage identification procedure 
then involves the full-field assessment of a suitable damage indicator. 

Effects of random error in the curvature data were also studied to ascertain the sensitivity of the method to 
noisy experimental data. The approach led to reasonably good predictions of the delamination shape and 
location even near the panel boundaries. 

Conclusions 

Two finite element based computational methods, SEA and iFEM, were reviewed and their salient 
features discussed in relation to SHM systems for use in future spacecraft, large space structures, and 
habitation structures. The methods and their associate algorithms satisfy the key requirements for use in 
real-time, large-scale structural applications, i.e., versatility, robustness, and computational efficiency. 
Several recent research efforts in which SEA and iFEM played a key role were also highlighted. These 
include computational and experimental studies of reconstructing, based on in-situ strain measurements, 
structural deformations (i.e., 'shape sensing’), identification of structural anomalies including 
delamitations in composite laminates, and sensor-position optimization studies employing genetic 
algorithms. The iFEM capability also permits reconstruction of internal structural loads - the essential 
information for mitigation of structural-failure accidents by means of cost- and time-effective service and 
necessary repair. 
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